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In this paper we attempt to reconstruct one of the last projects of Volodya Geyler 
which remained unfinished. We study motion of a quantum particle in the plane to 
which a halfline lead is attached assuming that the particle has spin ^ and the plane 
component of the Hamiltonian contains a spin-orbit interaction of either Rashba 
or Dresselhaus type. We construct the class of admissible Hamiltonians and derive 
an explicit expression for the Green function applying it to the scattering in such 
a system. 

1 Introduction 

There is no doubt that in the person of Volodya Geyler, science lost a bright 
personality, somebody who with safety and elegance treaded forward at the 
boundary between mathematics and physics, mastering the subtleties of the 
former and understanding deeply the meaning of the latter. For the authors 
of the present paper the sad news have a personal touch because the last 
talk he announced bore the title Exner-Seba hybrid plane with the Rashba 
Hamiltonian] he passed away at the opening of the conference in the Isaac 
Newton Institute in Cambridge where it had to be presented. 

We can only guess what Volodya intended to report but we decided that 
the best way to honour his memory is to reconstruct what we think could be 
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the contents of this lecture. We take the "hybrid plane" which we introduced 
twenty years ago and look what happens if the particle living in it has a spin 
and is subject to spin-orbit interaction in the plane part of the configuration 
space; we consider the two interaction forms which were objects of interest 
recently. We construct the class of admissible Hamiltonians described by 
the boundary condition in the coupling point between the plane and the 
halfiine lead attached to it. After that we derive an explicit expression for 
the respective Green's functions using standard Krein's function technique, 
and show how one can use them to derive properties of such system. 

Since this is not intended to be an in-depth study, we will speak mostly 
about the simple situation when there is no external field. It is easy to extend 
the results to the situation when the particle is subject to a homogeneous 
magnetic field in the plane. We will comment on this case briefly hoping that 
this and other ideas related to the problem will find a continuation. 

2 Spin-orbit interaction 

Let us first recall how one describes a two-dimensional particle with spin-orbit 
interaction; as announced we will pay most attention to the simple case where 
there is no external field, adopting the notation from |BGP07j . For a particle 
with two spin states the state Hilbert space is Tipiane = Iv^(R^, C^) and its free 
motion is described by Hq = 2^P^'^0; where pj = —ihdj, j = 1, 2, as usual, 
and (To is the 2x2 unit matrix. The spin-orbit interaction is introduced in 
different way: one is the so-called Rashba Hamiltonian 

i^R := ifo + ^f^R , ■= 0-1P2 - cr2Pi , (2.1) 

n 

where or G M is the Rashba constant and aj are the usual Pauli matrices, 
the other is the Dresselhaus Hamiltonian 

Hb := ^0 + ^f>D , f>D := a2P2 - aipi , (2.2) 

in which interaction strength is given by Dresselhaus constant od G M. 

Since the choice of the units will not be important in the following we 
get rid of the constants in the usual way introducing k := h~^p and Xj : = 
h~'^Tn*aj, j = R, D. Up to the multiplicative factor, Hj = -^^Hj, J = R, D, 
the both versions of the Hamiltonian acquire then the simple form 

Hj = Ho + 2xjUj , [/r := cxi/ca - o^afci , Ub := cTa/ca - aiki , (2.3) 
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with Hq := p^ctq, which we shall use in the following. 

As usual the properties of such a Hamiltonian are encoded in its resolvent. 
The latter is known explicitly from the paper |BGP07] . By a nice algebraic 
trick, so characteristic for the work of Volodya Geyler, the derivation is re- 
formulated as a scalar problem which involves the well know resolvent kernel 



G'o(x,x';2) = ^Ko( 



z|x — x'l) of the Laplacian in L^(R^), where Kq is 

jsion 

^ (2.4) 



the zero-order MacDonald function. It leads to the expression 

G'j(x,x'; z) 
with the diagonal elements 
G'f(x,x';;.) 



G'f(x,x';z) Gf(x,x';z) 



Gj^ ( X, X ; z 
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(^o(Cj^ 



^o(Cj- 



i^-{z + X2) 

x'|))+Ko(C/|x-x'|) + iro(Cj-|x-x'|) 



for both the J = R, D, while the off-diagonal ones are 



G]^(x, X z] 



1{X2 — Xo) — [Xi — X 



\z + x^) |x — X' 



u=± 



G^2(x,x';z) 

and Gj^(x, x'; z) 

as 



{X2 - X2) - l(Xi - Xi^ 



An Ia/— (^T^ijy |x — X' 



GYW,^;z 



the effective momenta are defined at that 



Cf:=^-iz + Hj)±i>cj (2.5) 

For the hybrid plane model which we will describe in the next section 
we will need also the renormalized Green's function, i.e. the diagonal value 
obtained after the subtraction of the divergent term. 



GTiz) 



lim 



Gj(x, x'; 2;) + — In |x - xVo 
Zir 



(2.6) 



notice that the limit is independent of the position x in view of the transla- 
tional invariance of the Hamiltonian Hj. By a straightforward computation 
|BGP07j one finds that the off-diagonal elements vanish in the limit while 



(^) 



(Q(c+)-Q(r)) + ^(Q(n + Q(r)) 



with Q{z) := 2^('0(1) ^ ^^^{~^) + 1^2) expressed through the digamma 
function, hence the renormahzed Green's function can be written as 



V.(l)-lln(-^)+ , + + 



(^0 ; 

(2.7) 

recall that ~ 0.577 in the above formula is Euler's number. 

The case when a homogeneous magnetic field B = perpendicular to 
the plane is applied is treated in a similar way. The momentum k in the 
Hamiltonian (12.31) has to be replaced with K = k — a where A = is the 
vector potential associated with the field, and the Zeeman term 76(73 with 
7 := hq^— has to be added. The the reduction to the scalar case works 
again and yields explicit expression for Green's functions |BGP07j in terms 
of confluent hypergeometric instead of Bessel functions. 

After this preliminary, we are ready to turn to our proper subject. 



3 Motion in the hybrid plane 

Let us now consider that the sy stem h as a mixed dimensionality and its 
configuration space consists, as in jES87j . of a plane described above to which 
a halfiine lead is attached; conventionally we place the junction to the origin 
of coordinates in the plane. As it carries the same spin | particle the lead 
component Hilbert space is Tiiead = L'^(M+,'C^), and the whole state space 
of the system is the consequently the orthogonal sum H := TC\ead © '^^piane- 
The wave functions are thus of the form \1/ = {■0iead, V'piane}'^ where each of 
the components is a 2 x 1 column. 

Our aim is to find the dynamics of the particle on the described configu- 



ration space. The idea of the construction is the same as in [ES87j . an d since 
it was used many times — let us just recall |BEG03l IBG031 IBGMP021 IES97 
as a sample of this work — we can just recall briefly the scheme and describe 
how it applies in the present situation. We start from the decoupled operator 
:= H\ead ffi Hj where the first component is the Laplacian on the halfiine 
i^iead'^iead = ""^lead with Ncumauu bouudary condition at the endpoint, and 
Hj is the Hamiltonian with the spin-orbit interaction discussed in the previ- 
ous section. We restrict to functions which vanish in the vicinity of the 
junction, obtaining thus a symmetric operator of deficiency indices (4, 4), and 
after that we seek admissible Hamiltonians among its self-adjoint extensions. 
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The problem differs from those mentioned above only by the presence 
the spin degree of freedom and the best way to characterize the extensions 
is again through boundary conditions. To this aim we need the boundary 
values. Those on the halfline are the columns '0iead(O+) and '0i'ead(O+)- 
plane the functions from the domain of the restriction have a logarithmic sin- 
gularity at the origin and the generalized boundary values i^j('0piane)5 j=0,l, 
appear as coefficients in the expansion 

^plane(x) = ^o(^plane) In |x| + Li(?/)piane) + o(|x|) . (3.1) 

ZTT 

Using this notation we can write the sought boundary conditions as 

Cad(0+) = ^^lead(0+) + C*Lo(^plane), 
^l(V'plane) = C-^lead (0+) + I^Lq (^/'plane) , 

where A, C, D are 2x2 matrices, the first and the third of them Hermitian, so 
the matrix A := ('^'^) depends of sixteen real parameters as expected. It is 
straightforward to check that the corresponding boundary form vanishes un- 
der the condition ( 13. 2p . and therefore each fixed A gives rise to a self-adjoint 
extension i/^ of the restricted operator. Notice that the analogous boundary 
conditions apply also to the magnetic case mentioned in the previous section 
due to the same character of the singularity. 

Let us further mention that the above boundary conditions are generic 
but do not describe all the extensions leaving out cases when the matrix A 
is singular; this fiaw can be corrected in the standard way - see, e.g., [AP05j 
- if one replaces (13. 2p by the symmetrized form of the relation, 

^/^lead(0+)\ ^^/^Ld(0 + ) 
V^o(V'plane)/ V-^l(^pl anc / 

where A, B are matrices such that (^|i3) has rank four and AB* is Hermitean. 
We will restrict ourselves, however, to the case B = —I in the following; the 
same is true for the alternative form of the b.c. mentioned below. 

The choice of the parameter matrix depends on the way in which the lead 
and the plain are connected. In particular, diagonal A, C, D correspond to the 
situation when the contact does not couple the spin states, and in addition, 
the matrices are scalar if the coupling is spin-independent. Moreover, the 
lead and the plane are decoupled if A is block-diagonal, i.e. C = 0. A naive 
interpretation is that C is responsible for the coupling while A and D are 
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point perturbations at the contact "from the two sides" . The reahty is more 
comphcated, though. If the halfline models a thin fibre of radius p coupled 
to the plane, then in the spin-indepedent case the natural choice seem to be 

^ = 7^ ctq , C = i — o-Q , D = -o-Q In p (3.4) 
2p \/2Tip 



in analogy with the discussion of then spinless case performed in |ES97 . 



4 The Green function 

Having defined the class of admissible Hamiltonians we can proceed to con- 
struction of their resolvents. One can use the standard procedure based 
Krein's formula - see, e.g., |AGHH] or |BGMPd2] . The starting point is 
Green function of the decoupled system which has a block-diagonal form, 

G(x,a;;x,x;.) = ^ Gj(x,x';.)J' ^^'^^ 

where G'j(x, x'; z) is given by ( 12. 4p and 

G'iead(a:, X] z) = cos ^/zx< Q-^^^> cTq 

with the conventional notation, x< := min{x,x'}, x> := max{a;,x'}, since 
we assumed Neumann boundary condition. The Krein function Q{z), which 
is an analytic 4 x 4-matrix valued function of the spectral parameter z, is 
defined through diagonal values of the kernel, with the above described renor- 
malization in the plane component, specifically 

- ( *o" GT(.) ) ■ <^-^' 

To express the full Green function it is useful to cast the conditions of the pre- 
vious section into an alternative form by changing the basis in the boundary 
value space: instead of the vectors appearing in (13.31) we consider 

-'plane , 



^o(V^plano) / V^l('0: 
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It is easy to see is that they satisfy Atii/j + 2'0 = with B = —I and 

V -CA-^ D - CA-^C* 

It is obvious that A = —AB* is Hermitean; we have to suppose, of course, 
that the matrix A is regular, or roughly speaking, that if_4 has no Neumann 
component on the halfline (notice that this is true, e.g., for f l3.4p ). The 
advantage of these boundary conditions is that our comparison operator 
is characterized by Tii/j = 0, i.e. A'^ = I, B^ = 0. In such a case we can 
use the result of |AP05] (which in our case boils down to the usual Krein's 
formula) by which the resolvent kernel of if_4 is given by 

G_a{x, X \ X, x'; z) = G^{x, x'; x, x'; z) (4.4) 

-G\x, 0; X, 0; z) [Q{z) - A] ~^G°(0, x'; 0, x'; z) 

differing from the free one by the second term on the right-hand side which is 
a rank sixteen operator. Notice that even if the coupling is spin-independent, 
A = ("^) ® do and similarly for A, the Green function does not decompose 
because spin states are coupled by the spin-orbit interaction in the plane. 




5 Properties of 

We will concentrate on the case when there is a nontrivial coupling between 
the two parts of the configuration manifold, i.e. A is no block-diagonal. In 
the opposite case we have two separate problems; the halfline one is trivial 
while spin-orbit Hamiltonians with point interactions deserve an investigation 
- we believe that the reader can find a study on this topic in the contribution 
of K. Pankrashkin to this issue. To keep things simple, we suppose that the 
coupling is spin- independent, ^ = (°^) ® do with c 7^ 0, so 

Let us first remark that the junction can bind. For instance, to any 
number in (— Xj,0) one can find i/_4 for which it is an eigenvalue. Indeed, 
writing the negative energy as — we see that (15. ip is singular if the relation 
(k"^ — a){G]'^'^{—K'^) — d) = |cp is valid, or in the original parameters 

(k - a) (G;^^'^(-k2) -d) = |c|2 . (5.2) 
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By (12.71) G'j^°(— K^) is real- valued for < xj, then it is easy to pick the 
parameters a, d in such a way that fl5.2p is satisfied. 

The most interesting aspect of the problem, of course, is the transport 
through the junction. A straightforward way to treat it is to use the formula 
(14.41) . Any vector of 7i can be written as [H^ — z)~^ijj^ for ijj^ G D[H^) and 
Imz 7^ 0, hence applying the formula to it we get 

= ^° - lAQi^) - A]-'^liH' - , (5.3) 

where jz '■ ^ H is the trace operator given by the kernel G'^(x, 0; x, 0; z) 
and 7* is its adjoint. Notice that T{z)*{H'^ — z)~^ip^ is just the vector of the 
values of at the connection point and Q{z)—A is position- independent, so 
the second term at the right-hand side is easy to compute. Now we employ 
the usual trick letting z to approach a real value /c^. The resulting function 
ceases to be L^, of course, but it still satisfies locally the boundary conditions 
in the junction and it can yield a generalized eigenfunction associated with 
the scattering which we are looking for. 

In particular, we can choose the vector with the "upper" component 
only, 'i/'piane = Aead = COS kx (uoticc that not every combination of 
gitifcx ^-jj since '^jead has to satisfy Neumann boundary condition at the 
origin). It is a straightforward exercise to invert the matrix fl5.ll) and to 
compute ; it yields the refiection amplitude of a particle travelling over the 
halfiine towards the junction with a momentum k in the form 

^ Iri - «) (^r(^^) - d) - |gp 



a) (Gf^'ik^) - d) 



cP 



naturally independent of the particle spin state. It is straightforward to 
recompute it in terms of the original parameters from ( 13. 2^ : we get 

{a + ik){Gr{k')-d) + \c\^ 
^ ' (a - ik){GT{k^) -d) + |c|2 ■ ^""^ 

Since Gf^{k'^) is generally complex |7^(A;)p is not equal to one for |c| 7^ 
which is natural because the coupling allows the particle to pass from the 
lead to the plane. In particular, in the absence of the spin-orbit coupling 
when the last term at the right-hand side o f (12.71) is missing, (15. 4p reduces 



to the refiection amplitude derived in |ES87j for spinless case (up to the sign 



of k which is due to the opposite halfiine orientation used there). 
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In the magnetic case one can proceed in the same way replacing G^f'^{k'^) 
by the renormahzed magnetic Green's function of [BGP07j . There is a sub- 
stantial difference, though. Now the Green function is real- valued, and conse- 
quently, the scattering on the halfiine is unitary, \Tl{k)\'^ = 1. The scattering 
in such a case will naturally exhibit resonances due to the discrete spectrum 
of the spin-orbit Hamiltonian in the plane which are worth of investigation, 
however, this is a subject for another paper. 
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